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I. INTRODUCTION 

In two recent pa^pers^^ hereafter denoted by I and II, 
respectively, an effective field theory has been developed 
that describes the quantum ferromagnetic transition in 
disordered itinerant electron systems, and the exact crit- 
ical behavior was determined for all spatial dimensions 
d > 2. Soft modes in addition to the order parame- 
ter fluctuations, viz., diffusive particle- hole excitations 
( "diffusons" ) , were found to play a crucial role in the de- 
termination of the critical behavior. This theory studied 
the transition as it is approached from the paramagnetic 
phase. The critical behavior in the ferromagnetic phase, 
in particular the exponent /3, was obtained from scaling 
arguments. This raises the following question: In the fer- 
romagnetic phase of a Heisenberg magnet, the existence 
of spin waves, which are the Goldstone modes that result 
from the spontaneous breaking of the rotational symme- 
try in spin space, changes the soft-mode spectrum com- 
pared to the paramagnetic phase. Given the importance 
of the soft modes for the critical behavior, will the Gold- 
stone modes invalidate the scaling arguments based on a 
theory for the paramagnetic phase? 

A related question follows from the fact that the Gold- 
stone modes, contrary to the soft particle-hole excita- 
tions, exist at nonzero temperature as well as at T = 0. 
Therefore, if there were an observable whose leading criti- 
cal behavior depended on the Goldstone modes, then one 
would expect that observable to exhibit classical critical 
behavior on the ferromagnetic side of the phase transi- 
tion, even at temperatures low enough for the system to 
otherwise exhibit quantum critical behavior. In such a 
case, the scaling arguments used in II might still be for- 
mally correct, but their region of validity would vanish. 

In the present paper we investigate the validity of the 
scaling arguments used before. We will find that they 
were correct, and the critical behavior obtained in II was 
indeed exact, even on the ferromagnetic side of the transi- 
tion. To this end, we develop the ferromagnetic analog of 



the theory described in I and II. The general strategy is to 
derive a theory that expresses the system in terms of the 
order parameter fluctuations, and any other soft modes 
that couple to them. We then analyze the resulting the- 
ory to determine the critical behavior of the system as it 
approaches the transition from the ferromagnetic phase. 

This paper is organized as follows. In Sec.^we present 
a simple generalized mean-field theory that allows an ex- 
plicit description of the ferromagnetic phase based on 
the theory developed in I. Its results agree with those 
obtained in II, except that it does not reproduce the loga- 
rithmic corrections to power-law scaling found in the lat- 
ter paper. In Sec. lIIII we present and motivate an effective 
action for itinerant quantum ferromagnets that describes 
the instability of the ferromagnetic phase. This action 
is analyzed by means of renormalization-group methods 
in Sec. IIVI It is shown that the generalized mean-field 
theory yields indeed the exact critical behavior except 
for logarithms. In Sec. [3 we summarize and discuss our 
results. Several technical points related to a derivation 
of the effective action from a microscopic model are rel- 
egated to an appendix. 



II. GENERALIZED MEAN-FIELD THEORY 

In this section we present a simple way to obtain an 
approximate theory for the ferromagnetic phase from the 
formalism of I. It will turn out that this simple theory 
produces essentially the correct result, and the techni- 
cally more involved development in the later sections will 
serve to show this. 



A. A simplified effective action 

We recall that the effective action AcS in I took the 
form of a Landau-Ginzburg- Wilson (LGW) functional 
of the order parameter field M, a generalized nonlinear 
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sigma-model for the soft fermionic degrees of freedom q, 
and a term that couples these fields, 



Aff = Agw[M] + Al^mM + Ac[M, q] 



(2.1) 



We simplify the expressions given in I for these terms 
by making a mean-field approximation for the order pa- 
rameter field, i.e., we replace the fluctuating vector field 
M by its expectation value, which is characterized by a 
number to, 
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Here (x, r) denote the dependence of M on real space 
and imaginary time, k and n are the wavevector and the 
Matsubara frequency index, respectively, of the Fourier 
transform of M, a is a replica index, V denotes the sys- 
tem volume, and we have assumed that the magnetiza- 
tion orders in the z or 3-direction. With this approxi- 
mation, the LGW part of the action becomes simply a 
Landau theory. 
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The coupling term in the action (Eq. (2.25c) in I) sim- 
plifies to 



Ac ^ m \f^t dxJ2Yl 

a r=0,3 
X^tr [{Tr(E>S3)Q^Ux) 



(2.4) 



Here Kt is the spin-triplet interaction amplitude of the 
underlying fermionic model; the spin-triplet interaction 
term has been decoupled by introducing the Hubbard- 
Stratonovich field M. Q is a composite field^ that origi- 
nates from the fermion variables^^ 



-V'U'02T 

\-i>i'\i>2'\ 



-V'iT'02i 



-V^l|V2i 

-V:uV'2i 



'01tV'2T \ 

-V^liV'2T 

-011^21 / 



(2.5a) 



Here the ip and ip £^re the fermionic, i.e., Grassmann- 
valued fields that provide the basic description of the 
electrons,^ and all fields are understood to be taken at 
position X. The indices 1,2, etc. denote the dependence 
of the Grassmann fields on fermionic Matsubara frequen- 
cies LOn-i = 27rr(ni -1-1/2) and replica indices ai, etc., and 
the arrows denote the spin projection. It is convenient to 
expand the 4x4 matrix in Eq. (|2.5al) in a spin-quaternion 
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(2.5b) 



with tq = sq = I2 the 2x2 unit matrix, and Tj = —Sj = 
—icTj, (j = 1, 2, 3), with (Ti_2,3 the Pauli matrices. In this 
basis, i — and i = 1,2,3 describe the spin-singlet and 
spin-triplet degrees of freedom, respectively. The r — 0,3 
components correspond to the particle- hole channel (i.e., 
products of i^ijj), while r = 1,2 describe the particle- 
particle channel (i.e., products of or V'V')- For our 
purposes the latter are not important, and we therefore 
drop the r = 1,2 from the spin-quaternion basis, as was 
done in I and II. In terms of the remaining fields, the spin 
density can be expressed as 

m ab 
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(2.5c) 



and the magnetization couples linearly to the spin den- 
sity, cf. Eq. H2.4|l . The difference between Q and Q is 
that the former contains all fermionic degrees of freedom, 
while the latter contains the soft fermion excitations only. 

The latter can be parametrized in terms of a matrix field 

„ 1.3 
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(2.5d) 



The elements of this block matrix, clockwise from up- 
per left, denote those elements of Q that have frequency 
labels ni,n2 > 0; rii > 0,n2 < 0; n.i,n2 < 0; and 
ni < 0,712 > 0, respectively. 

The nonlinear sigma-model can be expressed in terms 
of Q afsiiiS&i 
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+2H j dx tr (qQ{x)^ +A["J,[Q] 
with Q a frequency matrix with elements 

^12 = (to 8) So) S12 OJni ■ 



(2.6a) 



(2.6b) 



The coupling constants G and H are proportional to the 
inverse conductivity and the specific heat coefficient, re- 
spectively. A\^^ is the spin-singlet interaction. 
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, (2.6c) 



with a coupling constant Kg. At the level of the bare 
effective action, .Anlo-m does not contain any spin-triplet 
interaction term, since the Hubbard-Stratonovich decou- 
pling turned the latter into the terms containing M . For 
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later reference we note, however, that under renormaUza- 
tion a spin-triplet interaction will be generated as long 
as Kg 7^ 0. We therefore add to the action a term 



effects of the fermionic soft modes by formally integrating 
out q, 



A 



TTTKt 



d-Ei-^y E EE 

r— 0,3 ni,n2,m a i—1 



(2.6d) 



with Kt the spin-triplet interaction constant generated 
by renormalization. 



B. Generalized Landau theory 

We now can obtain a generalized mean-field or Landau 
free energy density, f{m), that takes into account the 
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with N the number of replicas. Since AcS contains q to 
all orders, the integral in Eq. H2.7I) can be performed only 
in terms of a loop expansion. The lowest-order term is 
obtained by expanding Q to second order in q. The inte- 
gration can then be done exactly. The Gaussian propaga- 
tors that are needed for this calculation have been given 
in I (the generalization to the case ^ 7^ is straightfor- 
ward, and can also be found in Ref . |a) , and we therefore 
only give the result, 
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By minimizing /(m) with respect to m, we obtain a gen- 
eralized mean-field equation of state. If we introduce 
suitable units, drop non-essential constants, and add an 
external magnetic field h, the latter can be written 

h — tm + u w? 

const ^ -A {k? + nn)^n 

(2.9a) 

with const > 0. An inspection of the integral shows that 
the leading term for small m is finite for d > 2, and 
that the leading nonanalytic m-dependence is given by 
j^((i-2)/2 generalized mean-field equation of state 

thus reads 

h^tm + vm'^/'^ + um^ , (2.9b) 

with u > 0. 

We see that for 2 < d < 6 the nonanalytic term dom- 
inates over the contribution. The exponents (3 and 
6, defined by m{h = 0) cx {~tY and m{t = 0) cx h^/\ 
respectively, for 2 < d < 6 in this approximation are 

(3^2/{d-2) , . (2.10) 

These values agree with the ones found in II, apart from 
logarithmic corrections to scaling that the generalized 
mean-field theory misses. Equation (|2.9a|) is also very 



similar to the effective equation of state that was ob- 
tained in Ref. 0, and it has the same qualitative proper- 
ties. Notice, however, that the derivation of Eq. (|2.9a() 
did not involve any divergent integrals, while the equiva- 
lent result in Ref. 1^ was obtained by resumming an infi- 
nite series of divergent terms. We will discuss the relation 
between these theories in more detail in Sec. 

While the above procedure provides a fast and sim- 
ple way to incorporate soft-mode effects into the equa- 
tion of state, and the agreement with prior results is 
encouraging, our mean-field approximation is of course 
uncontrolled. Furthermore, since it neglects the order 
parameter fluctuations, it does not shed any light on 
the role played by the Goldstone modes. In the follow- 
ing two sections we will therefore construct an effective 
field theory for the ferromagnetic phase, and perform a 
renormalization-group analysis to study the critical be- 
havior. We will find that the results of the generalized 
mean-field theory are indeed exact apart from logarith- 
mic corrections to scaling. 



III. EFFECTIVE FIELD THEORY FOR 
FERROMAGNETS 

In this section we present an effective field theory for 
the ferromagnetic phase in analogy to the one for the 
paramagnetic phase given in I. The structure of such a 
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theory can be deduced without a calculation by combin- 
ing various existing results for disordered itinerant ferro- 
magnets. In the interest of brevity and clarity, we pursue 
this route. We have also checked the result by means of a 
derivation from a microscopic model, using the methods 
of I and Ref . We summarize the salient points of this 
derivation in Ap pendix IXI a complete account has been 
given in Ref. HU 

In I it was shown how to derive an effective long- 
wavelength and low-frequency field theory for a disor- 
dered metal that approaches a ferromagnetic instability. 
We now give the corresponding action for a ferromagnetic 
metal, which we will motivate below. We denote the fluc- 
tuations of the longitudinal part of the order parameter 
field M about its expectation value m by SMf , and the 
transverse part by Mt. Similarly, we denote the spin- 
singlet and the longitudinal spin-triplet components of 
the fermionic field qhy qi, and the transverse spin-triplet 
components by qt- The action then has the general form 



nels are similar to those in a paramagnet. The spin- 
singlet and the longitudinal spin-triplet fermionic fields 
are known to be unaffected by the presence of a nonzero 
magnetization (see, e.g., Ref. 0); the vertex for qe thus 
has the same form as in I. The same is true for the longi- 
tudinal order parameter fluctuations. The distance from 
criticality, i, is actually magnetization dependent, but 
this is irrelevant for power counting purposes. The cou- 
pling between SM^ and qi is also unaffected. Notice that 
this coupling provides the leading frequency dependence 
of the SMi propagator, viz., a term proportional to il/k^ 
with fc a wavevector. We will refer to this propagator as 
the longitudinal magnon. That is, in Fourier space and 
for small wavenumbers and frequencies the longitudinal 
magnon propagator has the form 

{6MfXk)SMe{-k)) = 



\t\ + ai,d-2\k\d-^ + a,,2 fc2 + Gi cj ^ n/k^ 



Aff[M,9] = A'i[6Ah,qi\+Af[Mt,qt] 
+l^A[5Mi,,Mt,qi,qt] . 



(3.1a) 



Here the superscript G denotes the action at the Gaus- 
sian level, where the longitudinal and transverse fields 
decouple, separated into longitudinal and transverse con- 
tributions, represents contributions of higher than 
bilinear order in the fields. We now specify the various 
terms in this action in the same schematic notation is in 
I (cf. Eq. (3.8) of I), suppressing everything that is not 
necessary for power counting, and considering all fields 
as functions of real space position x and Matsubara fre- 
quencies. The Gaussian contributions arc then given by 
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(3.1b) 
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(3.1c) 



For this reason we have left the less leading direct fre- 
quency dependence (oc f2) out of the SMg vertex. 

The Eqs. 1)3. l|l and H3.2|l demonstrate an important 
point that was discussed in detail in I and II. Namely, 
there is more than one time scale in the problem, and 
hence more than one dynamical exponent z.^^ The dif- 
fusive dynamics of the fermions are described by an ex- 
ponent z = 2, while the critical z, which describes the 
critical dynamics of the order parameter field, has a dif- 
ferent value. This will be important in what follows. 

The soft-mode structure in the transverse spin-triplet 
channel is qualitatively different compared to the para- 
magnetic phase. The transverse fermion fields acquire a 
mass proportional to m, see, Ref. @ or Eq. (|2.4|l . The 
order parameter vertex, on the other hand, is massless 
due to the presence of ferromagnetic Goldstone modes. 
The structure of the fermionic vertex suggests that the 
magnetization can scale like a gradient squared. For the 
transverse magnetization vertex, we expect terms that 
represent the classic dispersion relatioufiS 



fl ~ mk 



(3.3a) 



as well as terms that refiect the nonanalytic magnetiza- 
tion dependence of the magnon effective massjii^ 



(3.3b) 



Here X is a generic coupling constant that is equal to 
Ks in the spin-singlet channel, Kt in the longitudinal 
spin-triplet channel, and in the transverse spin-triplet 
channels. We now motivate this form of the Gaussian 
effective action. 

First of all, the magnetic order breaks the rotational 
symmetry in spin space and forces us to distinguish be- 
tween the longitudinal and transverse components of all 
fields. The soft-mode structures in the longitudinal chan- 



Again the leading dynamics come from the coupling to 
the fermions, which here produces a term il/m due to the 
mass in the qt vertex. Eqs. 1)3. 3|l therefore show that we 
must include a simple gradient-sqared term, as well as a 
gradient-squared term with a coefficient proportional to 
^(<i-4)/2 -j-j^g vertex in order to accurately represent 
the Goldstone modes. The latter reflects the same non- 
analyticity as the d^~^ term in the longitudinal channel 
(recall that ni ^ d^). The transverse magnon or Gold- 
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stone mode thus has the form 
{Mt{k) Mt{-k)) = 

1 



(3.4) 



Beyond Gaussian order, the transverse and sin- 
glet/longitudinal spin channels are coupled. By analyz- 
ing the spin structure of the action as given in I, we find, 
schematically. 



C2,i%/T j dx 6Meqf +C2^2Vt J dx SMgq^ + C2.3VT J dx SMtqeQt 
+uiT J dx {5Mif + u2T j dx {SMe^ {6Mtf +U3T J dx {6Mtf 



dx 
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dx 
dx 
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dl + i/4,2 ^ + m 
■ dl + Hi.3 n + m 



4,3 



2 2 
<ll It 



(3.5) 



The first and the third class of terms are the Mq^ and 
vertices, respectively, that were found to be impor- 
tant in describing the critical behavior as the transition 
is approached from the paramagnetic phase. The Mq^- 
vertices would generate the nonanalyticities in the Af^- 
vertices if they had not been included in the bare action, 
and together with the g^-vertices they give rise to log- 
log-normal corrections to power laws. We also include an 
M^-vertex since its coupling constant, m, is a dangerously 
irrelevant operator with respect to the magnetization. 



IV. RENORMALIZATION-GROUP ANALYSIS 

We now conduct a power-counting analysis of the ac- 
tion given in Sec. IIIII The purpose of this exercise is to 
determine the role of the Goldstone modes in the ferro- 
magnetic transition, and thus compare the critical be- 
havior on the ferromagnetic side of the transition with 
the results from the paramagnetic theory. 

For the following analysis, we assign lengths L a scale 
dimension of [L] = —1. Rescaling lengths by a fac- 
tor h under a renormalization-group transformation will 
change all other quantities according io A 6^"*!^, 
with [A\ the scale dimension of A. In particular, the 
temperature T and frequency have a scale dimension 
[T] = [Vl] = z. The scale dimensions of the fields we 
characterize as usual by means of exponents 



[5M,] = (d-2 + 7,,)/2 
[Mt] = (d~2 + T^t)l'2 
[qi,t] = (d-2 + r;;.,)/2 



(4.1a) 
(4.1b) 
(4.1c) 



As in Sec. Illll we consider all fields functions of frequency 
and real space position. Throughout this section we use 
Ma's technique of choosing scale dimensions, and then 
checking self-consistently whether this choice leads to a 
physical fixed point ^ 



A. Stable fixed point 

Before considering the critical fixed point that de- 
scribes the ferromagnetic transition, it is illustrative to 
discuss the stable fixed point that corresponds to the sys- 
tem deep in the ferromagnetic phase. We start with the 
longitudinal degrees of freedom. 

The natural choice for t is to be marginal at the stable 
fixed point This choice determines the scale dimension 
of 5Mi, so we have 



M = o , 7?, = 2 



(4.2a) 



qi we expect to be unaffected by the magnetization, so we 
choose its scale dimension to be consistent with diffusive 
behavior. 



4 = 



(4.2b) 



By the same argument, the dynamics of q^, i.e., the scale 
dimensions of U, and T in the vertex, should be gov- 
erned by a dynamical exponent z = 2. We further choose 
the coupling q^i to be marginal, which implies that the 
dynamics of 5Mi, as represented by the factor \/T in the 
last term in Eq. Ij3.1bp . are also governed by a z = 2. We 
thus have 



M = 



2diff = Zl = 2 



(4.2c) 



From Eq. (|3.1bp we then find that Gt, He, and K are all 
marginal. 



[Ge] = [He] = [K] = 



(4.2d) 



while ae^d-2 and ae^2 are irrelevant with scale dimensions 
[ae^d-2] = —{d— 2) and [ae,2] = —2, respectively. 

We now turn to the transverse degrees of freedom, qt 
is expected to be a massive fluctuation deep in the ferro- 
magnetic phase, which means we choose 



(4.3a) 
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This, together with the marginality of K, imphes that 
TO is marginal, as one would expect on physical grounds, 
and so is h, since [h] = [to] + [t], see the equation of state, 
Eqs. (E3, 



[m] = [h]=Q 



(4.3b) 



while 1/Gt and Ht are irrelevant with scale dimensions 
[1/Gt] = [Ht] = —2. We further choose at^2 and ct^i to 
be marginal. This renders at,d-2 marginal as well. 



[at,2] = [at,d-2] = [ct,i] = 



and it leads to 



Vt = , 



(4.3c) 



(4.3d) 



and to a transverse time scale governed by a dynamical 
exponent 



(4.3e) 



We see that at the stable fixed point the various time 
scales coincide, and there is only one dynamical expo- 
nent, z = 2. This reflects the fact that, deep in the 
ferromagnetic phase, the diffusive electron dynamics, the 
frequency dependence of the longitudinal magnon prop- 
agator, Eq. (|3.2|l . and the dispersion of the Goldstone 
modes, Eqs. ()3.3|1 . all lead to the same scaling of the fre- 
quency with the wavenumber. From Eq. (|3.5|) it is then 
easily seen that all non-Gaussian terms in the action are 
irrelevant with respect to the stable fixed point, with the 
least irrelevant operators having scale dimensions equal 
to —{d — 2)/2. Temperature and frequency are relevant 
operators, of course. This establishes the description of 
the ferromagnetic phase within our formalism. 



B. Critical fixed point 

We now turn to the critical fixed point that describes 
the ferromagnet-to-paramagnet transition. Since this is 
a symmetry-restoring transition, any sensible candidate 
for the critical fixed point must have the feature that the 
longitudinal and transverse fields have the same scale di- 
mensions. We demand that the fermion fields be diffu- 
sive, 



Ve = ri't = , 
with a time scale given by 

^diff = 2 . 

Equations H3.1bl3.1c(l then imply 

[G,] = [Gt] - [H,] = [Ht] = [K] 







(4.4a) 



(4.4b) 



(4.5) 



We furthermore require that cu and ci^t are marginal, 
which implies 



[SMe] = [Mt] = l + id~ z)/2 



(4.6) 



Here z is the dynamical exponent associated with the 
factor of VT in either of the coupling terms in Eqs. 
H3.1bl3.1c(l . As was explained in I, this z can equal to 
ZdiS, e.g. in terms in perturbation theory where the lon- 
gitudunal magnon propagator is convoluted with a diffu- 
sive one, which makes the frequency in Eq. H3.2|l scale like 
a wavenumber squared, making the longitudinal magnon 
effectively massive. The critical longitudinal magnon, on 
the other hand, we expect to be massless. In this case, 
we choose a(_d~2 to be marginal for 2 < d < 4, and ai^2 
for d > 4. This implies a critical time scale 



_\d for 2 < d < 4 
I4 for d > 4 , 

as well as a critical exponent ij = rjg = r]t, 

_f4-d for2<d<4 
~ I for d > 4 . 



(4.7) 



(4.8) 



t is by definition the only relevant operator (apart from 
the temperature/frequency and the external magnetic 
field) at a physical critical fixed point, and its scale 
dimension determines the correlation length exponent 
= l/[t]. Wc thus have 



l/{d-2) ioi-2<d<4 
1/2 ford>4 . 



(4.9) 



We now turn to the scale dimension of to. From Eq. 
(|3.1c|l . in conjunction with Eqs. H4.4a|l and H4.5|l . we have 



[m] = 2 



(4.10a) 



However, in order to determine the physical order pa- 
rameter exponents P and S, we need to take into account 
that the critical behavior of the magnetization is affected 
by dangerous irrelevant variables. From the equation of 
state, Eq. lf2Jb|l . we see that m{h = 0) cx -y-^/C^'-s) 
for 2 < d < 6, and m{h = 0) oc u^^/^ for d > 64^ 
The effective scale dimension of to is therefore [to] off = 
[m] + 2[v]/{d - 2) for 2 < d < 6, and [TO]off = [m] + [u]/2 
for d > 6. But from the equation of state we have 
[u] ^ [t] - 2[m] = -2 and [v] = [t] - (d - 2)[to]/2. This 
yields [u] = -2 for d > 6, and [w] = for 2 < d < 4, and 
[ti] = 4 - d for 4 < d < 6. Therefore, 



TO eff 




2) 



for 2 < d < 4 
for 4 < d < 6 
for d > 6 , 



(4.10b) 



which leads to 

_j2/(d-2) for2<d<6 



\l/2 



for d > 6 



(4.10c) 
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Finally, the effective scale dimension of the external 
magnetic field h is 



[^]cff = M + [m]cs = 




which implies 



for 2 < d < 4 
for 4 < d < 6 
for d > 6 , 

(4.11a) 



d/2 for 2 < d < 6 
3 for d > 6 . 



(4.11b) 



An inspection of Eq. (|3.5|l shows that all corrections 
to the Gaussian action are irrelevant with respect to the 
Gaussian fixed point, except that the C2 are marginal in 
the event that the factor %/T in this coupling carries the 
diffusive time scale. It was shown in I and II that this 
can indeed happen, and that this makes C2 marginally 
relevant with respect to the Gaussian fixed point. The 
actual critical fixed point therefore contains the effects of 
C2. This is a result of the existence of two time scales in 
the problem. Furthermore, the terms oiO{q^), which are 
irrelevant by power counting, turn out to be effectively 
marginal as well. This was also shown in I and II, and the 
logarithmic corrections to scaling that result from these 
marginal operators were explicitly determined. 

There is no need to repeat the discussion of the loga- 
rithmic corrections to scaling, since it is now clear how 
this solution carries over to the present case. Above we 
have shown explicitly that scaling works in the ferromag- 
netic phase, with the magnetization having a scale dimen- 
sion [m] = 2. This means in particular that the results 
obtained in II from scaling arguments for the free energy 
were correct. For instance, the exponent (3 is given by 
the result of the generalized Landau theory, Eq. H2.10|l . 
with logarithmic corrections as given in Eq. (3.6e) of II. 



correct critical behavior in all dimensions d > 2 apart 
from logarithmic corrections to power-law scaling. The 
results of the Landau theory are the same as those orig- 
inally obtained from a nonlocal order parameter theory 
in Ref. (Oj 



B. Hertz's fixed point 

We briefly discuss how Hertz's fixed point™ relates to 
the above discussion. Suppose one ignored the mode- 
mode coupling effects that are represented by the coef- 
ficients and at.d-2 in the Gaussian action, Eqs. 
(|3.1|) , and by v in the equation of state, Eq. (|2.9b|) . Then 
one has, for all d > 2, 



[t] = [m] = 2 , [u] = -2 



which leads to 



[m]cff = 1 



[h]c 



(5.1a) 



(5.1b) 



The dynamical critical exponent is then 

^c = 4 , (5.2a) 

and all static exponents have mean- field values, 

r, = , = = 1/2 , 6^3 . (5.2b) 

Of course, for 2 < d < 4 this fixed point is unstable 
against the mode-mode coupling effects as was discussed 
in detail in I. For 4 < d < 6 it is actually stable, and 
the only reason why Hertz's theory does not yield the 
correct critical behavior is that it misses the leading dan- 
gerous irrelevant variable for the magnetization, which is 
V rather than u. For d > 6 the exact critical behavior is 
mean- field like. 



V. DISCUSSION AND CONCLUSION 

We conclude with a summary of our results, and some 
additional remarks. 



A. Summary 

In summary, we have constructed an effective theory 
for the instability of the ferromagnetic phase of a dis- 
ordered itinerant Heisenberg ferromagnet at the quan- 
tum critical point. We have shown that the presence 
of ferromagnetic Goldstonc modes, or spin waves, does 
not change the critical behavior compared to the one ob- 
tained by supplementing results from the paramagnetic 
phase with scaling arguments. We have also given a very 
simple generalized Landau theory for this problem, which 
takes into accout the effects of soft fermionic modes in- 
dependent of the order parameter, and which yields the 



C. General remarlcs 

We finally come back to some of the points mentioned 
in the Introduction. We have shown, by explicitly con- 
sidering the ordered phase, that the scaling arguments 
used in I and II to extract the critical behavior of the 
magnetization were correct, and can be justified by a RG 
analysis. In particular, the presence of magnetic Gold- 
stone modes does not invalidate these arguments since 
the transverse fermionic modes and the magnetic Gold- 
stone modes essentially just switch roles as one goes into 
the ferromagnetic phase. We mention, however, that we 
have not addressed the question of the size of the quan- 
tum scaling region for any particular observable. This 
question can only be answered by an explicit crossover 
calculation for the observable of interest, which follows 
the behavior of the observable through the entire critical 
region at T > 0. 

Another point of interest is the relation between the 
theory developed here and Ref. 0. Like Hertz's theory. 
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the latter was a pure Landau- Ginzburg- Wilson theory, 
i.e., a field theory formulated entirely in terms of the or- 
der parameter field; the fermionic degrees of freedom had 
been integrated out. This integrating out of soft modes 
resulted in nonlocal vertices that made the theory un- 
suitable for explicit calculations, but the critical behavior 
could be extracted by a combination of power counting 
and scaling arguments. It is interesting to see that the 
current series of papers has completely vindicated this 
treatment. In particular, the equation of state derived 
in Ref. 0was Eq. H2.9a|l expanded in powers of , and 
the critical behavior determined in the nonlocal theory 
was indeed exact except for the logarithmic corrections 
to scaling in 2 < d < 4 that were discussed in II. It is 
also interesting to see that the much simpler generalized 
Landau theory of Sec. IIIBI which is based on a Gaus- 
sian approximation, reproduces the results of Ref. 9, and 
thus is also exact except for the logarithms. The reason is 
that, as we have seen, only two classes of non-Gaussian 
terms ultimately contribute to the critical fixed point, 
namely, the terms of 0{Mq^) and those of 0{q'^), and 
both of these lead only to logarithmic corrections to the 
Gaussian critical behavior. Section ITU therefore provides 
a very simple way to obtain the essentially exact critical 
behavior, but of course this is clear only a posteriori once 
the RG analysis has been performed. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE ACTION 

In this appendix we give the exact effective action at 
the Gaussian level, as it emerges from a derivation from 
the microscopic theory. A schematic version of this re- 
sult that is sufficient for power counting was given in 
Eqs. H3.1|l . and motivated by general arguments in Sec. 
Unl We start with the general field theory, Eqs. (2.10) 
of I. The first step is to find a saddle-point solution that 
is appropriate for a ferromagnetic phase. One then sep- 
arates soft and massive modes, and integrates out the 
latter in tree approximation. The result is the desired 
effective field theory for the soft modes. 



1. Stoner saddle point 

Equations (2.10) of I allow for a homogeneous saddle- 
point solution that reproduces Stoner theory. This is the 
saddle point considered in Ref. Hol supplemented with a 



saddle-point value for the order parameter field M, 

iQl2{x)\sp = Si2[SroS,oQl,+Sr3Si3Ql^] , (Ala) 

^Ai2(a;)|sp = 5i2[-SroStoi^ni + SrsS.si Am] , 

(Alb) 

'Afi(a;)|,p = S,,,Q6^3m/VT . (Ale) 

Here we use the same notation as in I, with 1 = (ni, ai), 
etc., comprising both frequency and replica indices. By 
substituting Eqs. ljAl|l into the effective action A (Eqs. 
(2.10) in paper I), and using the saddle-point conditions 
6A/6Q = SA/6A = 5A/5M = 0, one obtains the saddle- 
point equations 

1/2 



TO 

Ql 

A„ 



TTiVFTcl ^ 

m 

TriVFTcl 

A -(2rt)^/^TO . 
1 



-)3 iuJrnO 
1 e 



with 



Here 



(A2a) 
(A2b) 
(A2c) 
(A2d) 
(A2e) 

(A2f) 

(A3a) 
(A3b) 

with ^fc — /2m — PL. T s^t are the interaction amplitudes 
proportional to K^ t that were defined in I. Upon substi- 
tuting Eqs. (|A2a|l and l|A2f|l into (|A2e|l . one recovers the 
saddle-point equations of Ref. 10, and hence Stoner the- 
ory. 

It is useful to define various transport and thermody- 
namic quantities in terms of these Green functions. We 
will need 



G°(fc) 
Glik) 



[QZ{k) + g-{k)] 
[g+{k)-g-{k)] 



are Green functions given in terms of 



^7 ^'^Jl\\i9n[k) + g^{k)) 



+ ^k'Qt{k)gi(k) 



(A4a) 



and 



1 



TrmV n 



lim > 



{gt[k) + gi[k)) 



— k^gt{k)gi{k) 

dm 



(A4b) 
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as well as 

k 

These quantities represent the Born approximation for 
various conductivities and densities of states in the split- 
band system of Stoner theory. They are generalizations 
of the analogous quanties defined in Ref.0. For a physical 
interpretation of these quantities, see Ref. ITU 

2. Gaussion soft-mode theory 

The separation of soft and massive modes works in 
analogy to I, although the procedure is more cumbersome 



in the presence of a nonzero magnetization. The massive 
modes are integrated out in tree approximation to arrive 
at an effective soft-mode action, and the soft fermion 
modes are expanded in powers of q (cf. Eq. (|2.5dp ). again 
in analogy to I. The complete procedure and result can 
be found in Ref. 0. Here we list the Gaussian (i.e., 
bilinear in M and q) contribution to the effective action. 
The higher order terms are obtained analogously. 
The Gaussian action has the general form 

A^[q,M] = A^i^aM[q] + Ai^Gw[M] + Ac[q,M\ . (A5) 

We expand the q matrices in a spin-quaternion basis, see 
Eq. (|2.5b|) . The fermionic part of the action, which is the 
Gaussian part of a generalized nonlinear sigma model, is 
given by 



A 



NLcrM 



M - -4E E E E E 

a/3 ni,n2,n3,n4 k r.s— 0,3 



-5n, - 



Sa0 TTsTT' 



Here the n represent traces in spin-quaternion space. 



.0 
t 



= 2*r(r,.Tt)tr(s^s]) , = ^tr (raTrrj) tr (sgSis]) -f tr (s^sasj) 



--tr (Trr^^)tr {s^s^s^s]) 



(A6) 

(A7a) 
(A7b) 



The coupling constants G and H in Eq. ljA6p . as well as the densities of states iVp, are magnetization dependent 
generalizations of the analogous quantities in I. They are given by 



m 



16AK++ao-) , G^ ^lQ/n{a+ -a^) , G' ^ l&lTi{a+ + a^) , = 16/^(^0+ - a^-) ,(A8a) 



(iV+ + iVp)/2 , 7V^ = (7V+-7Vp)/2 



(A8b) 
(A8c) 



As can be seen by a direct comparison, the schematic action given by Eq. H3.1|l reflects all qualitative features of the 
complete result. Notice that the mass, i?A, in the transverse vertex comes from the coupling between M and q, cf. 
Eq. I|2.4(l . but has been included in y^NLo-M since it contains only fermionic fluctuations. 
The LGW part of the action that results from integrating out the massive modes is 



^LGWM - -EE(2-'5„o)E E Y^^+ {S.,^5,,+5.a5,2-i^^\2,^J)TtTY^{£t.+n,„r + {~f^'£;n+n^^^ 
a n>0 k = l \ m 

1 2 ■ 

'(fc)W«„(-fc). 

(A9) 



Here, we have adopted the notation of Ref. 10; V and nal and transverse spin channels, respectively. They are 
£ correspond to modified "diffusons" in the longitudi- 
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given by 



T^trnik) - 0±„(fe)/[l-0±„(fc)/27r7VFr] ,(AlOa) 
^Lik) = ?7±„(fc)/[l-77±™(fc)/2^^^Fr] ,(A10b) 



with 



and 



iVnm + Vira) ± {<ra + Vnra) , (AlOc) 
{vZ^-vfra)±{vl^a-vZa) ,(A10d) 



(AlOe) 

with G^'^ given by Eqs. ljA3|l . We have also introduced 
the notation 



T Q{m{m + n)) 



(All) 



The LGW vertex simphfies in the hniit of long wave- 
lengths and small frequencies. By usingi^i 



T5]l?±„(fe = 0) 



and 



T^f±™(fc = o) = -i/2r, , 



(A12a) 



(A12b) 



one recovers the functional form given in Eqs. (|3.1|l . ex- 
cept for the nonanalytic terms that emerge only at one- 
loop order and were added in Sec. IIIII 

Finally, the coupling between M and q to bilinear order 

is 

3 3 



12 k r,s=0,3 i=l j=0 

(Al3a) 

Here, 6 is a field with components 

n 

+ {-r+^^A'ri{k) . (A13b) 



Again, this has the same functional form as the schematic 
representation in Eqs. H3.1(l . 



3. Gaussian propagators 

The Gaussian propagators can be obtained by invert- 
ing the quadratic form given by the Gaussian action in 
the previous subsection. We do not list the complete 
propagators here, which can be found in Ref. Tl', but 
give only the hydrodynamic parts of the diffusons, which 
are needed in the calculation. For jim < one finds 



2itN, 



± 



(Al4a) 



Here, = 1/G±i?± with 



1/G^ 



1/G" ± 1/G^ 



(A14b) 
(Al4c) 



can be interpreted as the Boltzmann diffusivities in the 
upper and lower Stoner band, respectively. 

Similarly, one finds for rt > and m < 0, 



27riVF 



(Al5a) 



Here, = D/{1 ± 2iAr) = l/G±iJ, with 

1/G± = 1/G'± 1/G^ . (A15b) 



Notice that Eq. (|A15al) holds only for n > and m < 0. 
The result for n < 0, m > is given by the relationi^ 



(Al5c) 



By means of these expressions, one readily finds that 
the propagators have functional forms as given in Sec. IIIII 
Again, the nonanalytic dependences on k and A are not 
included in the Gaussian theory, and have been added in 

secirm 
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We use the notation, a ~ b for "a approximately equals 6" , 
a = b for "a is isomorphic to &", and a ~ 6 for "a scales 
like 6". 

z is defined by the scaling relation between frequency and 
wavenumber, Q ~ |fe|^. 

We note that it is permissible at this point to use the gener- 
alized mean-field equation of state, whose validity we want 
to establish. The point is that the equation of state serves 
only to obtain the values of the exponents /3 and S from [m], 
while the fact that [m] = 2 suffices to justisfy the scaling 
arguments used in I, which in turn establishes the validity 
of the equation of state up to logarithmic corrections. See 
also the last paragraph of this section. 
This distinction between the cases n > 0, m < and n < 0, 
m > was not made in Ref. 



